Local properties of local multiplicity distributions in hadronic Z decay * 
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Preliminary results on local multiplicity fluctuations in hadronic Z decays are presented. The data were obtained 
using the L3 detector at LEP. It is investigated to what extent Monte-Carlo models, which are tuned to reproduce 
global event-shape variables and single-particle inclusive distributions, can describe the local fluctuations measured 
by means of bunching parameters. 



1. INTRODUCTION 

The investigation of the evolution of the prob- 
abilities Pn{S) of detecting n particles in ever 
smaller sizes S of phase-space windows (bins) pro- 
vides detailed information on QCD multihadron 
production beyond single-particle densities, with- 
out trivial constraints from charge- and energy- 
momentum conservation. A deviation of this dis- 
tribution from that expected for purely indepen- 
dent particle production can be attributed to dy- 
namical local multiplicity fluctuations. 

The important quest behind such a study is the 
understanding of the origin of short-range corre- 
lations between final-state particles. As a conse- 
quence of these correlations, the normalized fac- 
torial moments (NFMs) Fq{S) of the local mul- 
tiplicity distribution exhibit a power- like 
increase with decreasing 6, namely Fq{S) oc 

, where 0^ are constants called intermittency 
indices. This phenomenon reflects the peculiar- 
ity of Pn{S) to become broader with decreas- 
ing S. Since NFMs satisfy the scaling prop- 
erty Fq{Xd) = X^'^'>Fq(6), this is widely regarded 
as evidence that the correlations exhibit a self- 
similar underlying dynamics. 

Local fluctuations in e+e~-processes have al- 
ready been studied is several experiments 1^ . The 
data do exhibit approximate power-like rise of the 
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NFMs with a saturation at small 5. The conclu- 
sion has been reached that such a phenomenon is 
a consequence of the multi-jet structure of events, 
i.e., groups of particles with similar angles result- 
ing in spikes of particles as seen in selected phase- 
space projections. Parton showers, fragmenta- 
tion, resonance decays and Bose-Einstein inter- 
ference can all contribute to these correlations. It 
has been found that for the statistics used at that 
time current Monte-Carlo models can, in general, 
describe the data, even without additional tun- 
ing. 

Recently, it has been realized that the factorial- 
moment method poorly reflects the information 
content of local fluctuations, since the NFM of or- 
der q contains a trivial contamination from lower- 
order correlation functions (see reviews |^). As a 
result, rather different event samples can exhibit 
a very similar behaviour of the NFMs. The fact 
that subtle details in the behaviour of Pn{S) are 
missing, together with the small statistics used, 
may be the reason why different Monte-Carlo 
models can reasonably describe the local fluctua- 
tions measured in e+e" annihilation so far. 

Another shortcoming of the factorial-moment 
measurement is that in moving to ever smaller 
phase-space bins, the statistical bias due to a fl- 
nite event sample (A'cv 7^ oo) becomes significant, 
especially for high-order moments q. This is be- 
cause in actual measurements the NFMs at small 
bin size are determined by the first few terms in 
the definition of the NFMs. In most cases this 
leads to a significant underestimate of the mea- 







Table 1 

NFMs and BPs for the distributions quoted. 



Distribution P„ NFMs BPs 

Positive-binomial C^p"(l-p)^-« nLi(l-^) 

Poisson exp(— p)/ri! 1 1 

Negative-binomial £0i^pn^p + i)-in+k^ UUi^ + i) 

Geometric nLi(l + «) ^ 



sured NFMs with respect to their true values. 

Cumulants are a more sensitive statistical tool 
(see and references therein). However, their 
measurement is rather difficult and was rarely at- 
tempted. Besides, the cumulants are expected to 
be influenced by the statistical bias to even larger 
degree, since they are constructed from the fac- 
torial moments of different orders q. 

2. LOCAL PROPERTIES 

An important step towards an improvement of 
experimental measurements of the local multiplic- 
ity distribution was made in where it was 
shown that any complex distribution can be rep- 
resented as 



n — -i+l 



(-5), 



i=2 



where A = Pi{S)/Po{S). The factor L„ measures 
a deviation of the distribution from a Poisson, 
for which Ln — 1. Non-poissonian fluctuations, 
therefore, exhibit themselves as a deviation of i„ 
from unity. The L„ are constructed from the 
bunching parameters (BPs) 



q Pq{S)P,.2{6) 



q>l. 



(1) 



The expressions for the BPs and NFMs for some 
popular distributions are shown in Table ^ The 
most interesting observation is that while the 
NFM is an "integral" characteristic of the P„ {6) 
and the BP is a "differential" one, both tools 
have values larger than unity if the distribution is 
broader than a Poisson. Generally, however, one 



should not expect that all BPs are larger than 
unity for a broad distribution; BPs probe the dis- 
tribution locally, i.e. they are simply determined 
by the second-order derivative of the logarithm of 
Pn{S) with respect to n. Note that in the case of 
local distributions, the width of the distribution 
is mainly determined by ri2{5). This observation 
is based on the simple fact that Pn{S) ceases to 
be bell-shaped at sufficiently small S. 

BPs are more sensitive to the variation in the 
shape of Pn{S) with decreasing S than are NFMs 
[^. In the case of intermittent fluctuations, one 
should expect 772 (^) oc (5"'*^. For multifractal lo- 
cal fluctuations, the riq{S) are (5-dependent func- 
tions for all g > 3, while for monofractal be- 
haviour rjq (S) = const for g > 3 Q . 

From an experimental point of view, the BPs 
have the following important advantages ^ : 

1) They are less severely affected by the bias 
from finite statistics than the NFMs, since the 
qth-order BP resolves only the behaviour of the 
multiplicity distribution near multiplicity n — q— 

1; 

2) For the calculation of the BP of order q, one 
needs to know only the g-particle resolution of the 
detector, not any higher-order resolution. 

In this paper, we present an experimental in- 
vestigation of local fluctuations in the final-state 
hadron system produced in decays at — 
91.2 GeV. The final-state charged hadrons have 
been recorded with the L3 detector during the 
1994 LEP running period. The calculations are 
based on approximately l.OM selected hadronic 
events. We compare the data with the JET- 
SET 7.4 PS §, ARIADNE 4.08 || and HERWIG 



5.9 ||9| models. The models have been tuned to 
the L3 data (lO). 

3. ANALYSIS 

1) Horizontal BPs: 

In order to reduce the statistical error on the 
observed local quantities when analyzing experi- 
mental data, we use the bin-averaged BPs : 



q- 



i(M) 



(2) 



where 7V,(M) = ijEra=iN,{m,S), N,{m,6) = 
being the number of events having q particles 
in bin m and M = A/J is the total number of 
bins (A represents the size of full phase-space vol- 
ume). To be able to study non-flat distributions, 
as for rapidity, we have to carry out a transforma- 
tion from the original phase-space variable to one 
in which the underlying density is approximately 
uniform, as suggested by Bialas, Gadzinski and 
Ochs 0. 

2) Generalized integral BPs: 

To study the distribution for spikes, we will 
consider the generalized integral BPs ^ using 
the squared pairwise four-momentum difference 
Q^2 = ~{Pi — P2Y ■ In this variable, the defini- 
tion of the BPs is given by 



q n,(Q2)n,_2(Q2) 



IP 

q- 



(3) 



where Ilq{Q^) represents the number of events 
having q spikes of size in the phase-space of 
variable Q12 1 irrespective of how many particles 
are inside each spike. 

To define the spike size, we shall use the so- 
called Grassberger-Hentschel-Procaccia counting 
topology for which a many-particle hyper-tube is 
assigned a size corresponding to the maximum 
of all pairwise distances (see Fig. |^). For a Pois- 
sonian production, the BPs are equal to unity 
for all q. 

3.1. In the rapidity variable 

In order to study fluctuations inside jets, in 
most investigations the fluctuations have been 
measured in the rapidity y defined with respect 




Figure 1. Example for the use of the Grassberger- 
Hentschel-Procaccia counting topology in two di- 
mensions. Full points illustrate the position of 
particles in an event. Points whose interparti- 
cle distance is smaller than or equal to the maxi- 
mum Q^, are connected by lines. This configura- 
tion corresponds to four spikes: one four-particle 
(1,2,3,4), two three-particle (3,4,5), (6,7,8) and 
one two-particle (9,10) spikes. 

to the thrust or sphericity axis (|. The analy- 
sis for this variable is performed in the full ra- 
pidity range | 1" |< 5. Fig. ^ shows the results 
for the BPs (H) for the rapidity variable after 
the Bialas-Gazdzicki-Ochs transformation. The 
second-order BP decreases with increasing M up 
to M ~ 20, which is found to correspond to the 
value of M at which the maximum of the mul- 
tiplicity distribution P„((5) first occurs at n — 0. 
At large M, all BPs show a power-law increase 
with increasing M, rjq ~ M"'. This indicates 
that the fluctuations in y defined with respect to 
the thrust axis are multifractal scale invariant. 
Note that the conclusion that fluctuations have 
a multifractal structure is possible without the 
necessity of calculating the intermittency indices 
4)q. In contrast, to reveal multifractality with the 
help of the NFMs, one first needs to carry out fits 
of the NFMs by a power law. 

Both JETSET and ARIADNE agree well with 
the data for the higher orders. However, ARI- 
ADNE underestimates 772 for all bin sizes while 
JETSET is significantly too low only for wide 
bins. HERWIG predictions (dashed lines) sig- 
nificantly overestimate the second-order BP ob- 
tained from the data. Since the second-order BP 
is determined by the dispersion of the distribu- 
tion, this means that the HERWIG produces too 
broad local multiplicity distributions. Such a re- 
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Figure 2. BPs as a function of the number M 
of bins in rapidity defined with respect to the 
thrust axis. The shaded areas represent the sta- 
tistical and systematic errors on the JETSET pre- 
dictions. 
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Figure 3. The second-order BP as a func- 
tion of the number M of bins in rapidity defined 
with respect to the thrust axis for hke-charged 
and unhke-charged particle combinations. The 
shaded areas represent the errors on the JETSET 
predictions. 



suit confirms that obtained by the ALEPH Col- 
laboration [Q. 

To study the disagreement in more detail, we 
split r/2 into two BPs: 



V2 



(±±) 

V2 



(4) 



Here, ?72^*'' is defined by (||) with N2{m,S) — 
N2^^\m, (5y), where N2^^\m, Sy) is the number 
of events having like-charged two-particle com- 
binations inside bin m of size 6y. Analogously, 
?72^ •* is constructed from the number of events 
iV2^ \m,5y) having unlike-charged two-particle 
combinations. Note that due to a combinatorial 
reason, ?72^^'' < ^ 

Fig. H shows that rf:^^'' and 772^ indeed be- 
have completely differently. While 772^^'' has the 
expected rise, 772"*^ ^ shows a strong decrease at 
low M and an onset of increase only at large M . 
The structure of 772 observed in Fig. ^ is a com- 
bination of these two effects. 

Let us remind that, in order to model the BE 
interference in JETSET, the momenta of identical 
final-state particles are shifted to reproduce the 
expected two-particle correlation function. The 
main disadvantage of such an ad hoc method is 
that it spoils the overall energy-momentum con- 
servation thus makes it necessary to modify also 
the momenta of non-identical particles to com- 
pensate for this. This effect can be seen in Fig. ||: 
JETSET without the BE correlations agrees with 
unlike-charged particles. However, JETSET with 
the BE effect disagrees with the data for unlike- 
charged particle combinations, while it agrees 
with like-charged. 

The strong anti-bunching tendency seen for 
unlike-charged particles at M < 30 can be at- 
tributed to resonance decays and to chain-like 
particle production along the thrust axis, as ex- 
pected from the QCD-string model . The lat- 
ter effect leads to local charge conservation with 
an alternating charge structure. Evidence for 
this effect was recently observed by DELPHI |Q . 
As a result, there is a smaller rapidity separa- 
tion between unlike-charged particles than be- 
tween like-charged and ^ is much larger than 



2^*'' at small M. Having correlation lengths 



(5y ^ 0.5 — 1.0 in rapidity, the resonance and the 
charge-ordering effects, however, become smaher 
with increasing M . 

3.2. In the four-momentum difference 

The study of BPs described above can help us 
to understand a tendency of the particles to be 
grouped into spikes inside small phase-space in- 
tervals. Another question is how the multiplicity 
of these spikes is distributed from event to event 
if the spike size goes to zero. To study this, we 
will use the BPs defined in (|). 

Fig. 1^ shows the behaviour of Xq ^ ^ func- 
tion of — InQ^. The full lines represent the be- 
haviour of the BPs in the Poissonian case. In 
contrast, all BPs obtained from the Monte Carlo 
models rise with increasing — In (decreasing 
Q^). This corresponds to a strong bunching ef- 
fect of all orders, as expected for multifractal fluc- 
tuations. The anti-bunching effect {xq < 1) for 
small — InQ^ is caused by the energy-momentum 
conservation constraint ||^. 

To learn more about the mechanism of mul- 
tiparticle fluctuations in the variable Q12, we 
present in Fig. |^ the behaviour of the second- 
order BP as a function of — InQ^ for multiparti- 
cle hyper-tubes (spikes) made of like-charged and 
those of unlike-charged particles, separately. For 
this kind of observables, HERWIG gives better 
agreement than the LUND models. 

4. DISCUSSION 

Local multiplicity fluctuations were studied by 
means of bunching parameters. Since all high- 
order BPs show a power-like rise with decreasing 
size of phase-space interval, none of the conven- 
tional multiplicity distributions given in Table | 
can describe the observed local fluctuations. 

For e+e~ interactions, one can be confident 
that, at least on the parton level of this reac- 
tion, perturbative QCD can give a hint for the 
understanding of the problem. Analytical calcu- 
lations based on the DLLA of perturbative QCD 
show that the multiplicity distribution of par- 
tons in ever smaller opening angles is inherently 
multifractal (l5j. Qualitatively, this is consistent 
with our results on the BPs for rapidity. Quan- 
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Figure 4. Generalized integral BPs as a function 
of the squared four-momentum difference be- 
tween two charged particles. The shaded areas 
represent the errors on the JETSET predictions. 
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Figure 5. Generalized second-order BP as a func- 
tion of the squared four-momentum difference 
between two charged particles. The shaded areas 
represent the errors on the JETSET predictions. 



titatively, however, the QCD predictions disagree 
with the e+e" data and MC models [l^ . 

In this paper we show that the power-law be- 
haviour of BPs is mainly due to like-charged par- 
ticles. JETSET gives the same power-law trend 
even without the BE effect. This means that the 
intermittency observed for like-charged particles 
appears to be largely a consequence of QCD par- 
ton showers and hadronization. 

The predictions of the ARIADNE 4.08 model 
are comparable with those of the JETSET 7.4 
PS model. This is essentially due to the same 
implementation of hadronization, which is based 
for both models on string fragmentation. 

A noticeable disagreement is found between 
the data and HERWIG model for rapidity vari- 
able. The conversion of the partons into hadrons 
in LUND models is based on the Lund String 
Model However, the hadronization in HER- 
WIG is modelled with a cluster mechanism p. 
This difference can be a rather natural candidate 
to explain the observed disagreement. A par- 
ticular concern is the large discrepancy for 772- 
The behaviour of 772 for not very small intervals 
is sensitive to low-multiplicity events, for which 
hadronization details could play a significant role. 
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